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We present preliminary results on the determination of spin-dependent, or polarised, Par- 
ton Distribution Functions (PDFs) from all relevant inclusive polarised DIS data. The 
analysis is performed within the NNPDF approach, which provides a faithful and statisti- 
cally sound representation of PDFs and their uncertainties. We describe how the NNPDF 
methodology has been extended to the polarised case, and compare our results with other 
recent polarised parton sets. We show that polarised PDF uncertainties can be sizeably 
underestimated in standard determinations, most notably for the gluon. 

The interest in spin-dependent, or polarised, Parton Distribution Functions (PDFs) of the 
nucleon is mainly motivated by the desire to understand its spin structure in terms of its 
quark and gluon parton substructure. It largely originates from the first EMC results pQ, 
originally interpreted as an indication that quark and anti-quark intrinsic angular momenta 
only contribute a small fraction of the full nucleon spin. A faithful knowledge of polarised 
PDFs is also an essential ingredient for exploring QCD beyond the helicity-averaged case and 
for studying the phenomenology of spin-dependent processes. 

Polarised PDFs have been investigated with increasing precision in recent years. On one 
hand, several experiments have contributed a large amount of data for a variety of processes, 
mainly inclusive polarised deep-inelastic scattering (DIS) but also proton-proton collisions and 
semi-inclusive reactions. On the theoretical side, the interest has been ultimately focussed 
on the global reconstruction of PDFs, together with their uncertainties. At least four groups 
have constructed such polarised PDF sets recently: BB [2], AAC [3], LSS [4] and DSSV 0. 
These sets slightly differ in the choice of datasets, PDF parametrisation and details of the QCD 
analysis (such as the treatment of higher-twist corrections). Nevertheless, they are all based on 
simple functional forms of the momentum fraction dependence of the PDFs at the refence scale 
(typically, a power-like behavior is assumed both at large and small momentum fraction) and on 
the Hessian approach for the estimate of uncertainties. Two main shortcomings are known to 
affect this methodology. The first one concerns how to propagate errors consistently from data 
to fitted parameters and then to observables: this is usually done by assuming Gaussian linear 
error propagation, which is not always adequate, in particular in those kinematical regions where 
few data are available. The second one consists in assessing the theoretical bias introduced by 
a fixed functional parton parametrisation. This is particularly delicate for polarised PDFs, 
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owing to the quantity and the quality of the data, which are respectively less abundant and less 
accurate than their unpolarised counterparts. 

In order to overcome these difficulties, in recent years the NNPDF collaboration has devel- 
oped a new approach to parton fitting (see, for example, [7J |H1 M HO] and references therein) . 
This new technique, designed to provide a faithful representation of PDFs and their uncer- 
tainties, is based on robust set of statistical tools, including Monte Carlo and Neural Network 
methods. 

In the NNPDF approach, experimental data are sampled by generating an ensemble of Monte 
Carlo replicas with data probability distribution; individual replicas are allowed to fluctuate 
in such a way that the mean value, standard deviation and correlation computed over Monte 
Carlo ensemble reproduce the experimental values, provided the sample is sufficiently large. 
Fitting an ensemble of parton distributions automatically propagates statistical fluctuations to 
the PDFs and then to observables. Hence, expectation values and uncertainties of PDFs (or of 
any observable) are obtained by considering their Monte Carlo integrals over the ensemble of 
replicas. Furthermore, in this approach neural networks are used as unbiased interpolants for 
PDF parametrisation. Since they provide functions depending on a large number of parameters, 
they are very flexible tools: this flexibility allows one to reduce the bias associated to the choice 
of some fixed functional form. 

The NNPDF approach has been succesfully applied to the determination of unpolarised 
PDFs and these NNPDF sets are routinely used by Tevatron and LHC collaborations for data 
analysis and data-theory comparisons. We will present here some preliminary results obtained 
by extending the NNPDF approach to the determination of a set of polarised PDFs. After 
illustrating the main features of our analysis, we will compare our results to those obtained 
by other collaborations. Specifically, we will see that the uncertainty on some polarised PDFs, 
most notably on the gluon PDF, are rather larger than previously estimated. 

The first NNPDF analysis of polarised PDFs, NNPDFpoll .0 henceforth, is based on a compre- 
hensive set of polarised DIS data. We exclude from our analysis data points with Q 2 < Q 2 ut — 
1 GeV 2 , since below such energy scale perturbative QCD cannot be considered reliable. We also 
impose W 2 > W 2 ut = 6.25 GeV 2 for the squared invariant mass W 2 = Q 2 (l — x)/x, according 
to the study presented in Ref. [TT]. This choice removes the dependence of results on possible 
dynamical higher-twist effects, which we do not include even though we do include target-mass 
corrections. The dataset used in the NNPDFpoll.O analysis is shown, after kinematic cuts, in 

Fig -B 

The experimental data used in this fit do not allow a full separation of individual flavour 
and anti-flavour parton densities. Hence, we define, for each light flavour q, the net amount of 
quark-antiquark spin density 

Aq(x, Q 2 ) = q^(x, Q 2 ) + gttfo Q 2 ) - q^[x, Q 2 ) + q^(x, Q 2 ) , 

where the superscript (tl) denotes that the parton spin is parallel (antiparallel) to the 
proton spin. We parametrise PDFs at the scale Qq — 1 GeV 2 by choosing, besides the gluon 
density Ag(x, Qq) = — g^, the following three linear combinations of light quarks: the 
flavour-singlet 

A2(x, Ql) = Au(x, Ql) + Ad(x, Q 2 ) + As(x, Q 2 ) , 
the non-singlet triplet and the non-singlet octet 

AT 3 (x,Q 2 ) = Au(x,Q 2 )-Ad(x,Ql) , AT a (x,Q 2 ) = Au(x,Q 2 ) + Ad(x,Q 2 ) -2As(x 1 Q 2 ) . 
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Each of these four combinations is parametrised by a neural network, with a total number of 
0(200) parameters, to be compared to O(10 — 20) used in other existing fits. 

A fast and accurate evaluation of polarised parton distributions, as required by the fitting, 
is achieved with the FastKernel method [TJ]. The accuracy of polarised PDF evolution has 
been shown to be O(10~ 5 ) comparing with the HOPPET code. 

Theoretical constraints are taken into account during the fitting procedure. We have im- 
posed positivity of physical cross-sections, which implies that the polarised structure function 
gi is bounded by its unpolarised counterparts Fx, so that \gi(x, Q 2 )\ < Fi(x,Q 2 ) [IS . For con- 
sistency, the unpolarised structure functions have been computed from the recent NNPDF2.1 
unpolarised PDF determination [10]. We have also used SU(3) symmetry to relate the first 
moments 

a 3 = dxAT 3 (x, Qo) > a 8 = / dxAT 8 (x, Ql) , (1) 

Jo Jo 

to the determination of et 3 and a g from baryon decay constants (allowing for large uncertainties). 
We have also performed a fit in which we 
have relaxed the first constraint in Eq. [T] and 
we have considered a 3 as a fit parameter, in 
which case we have found a 3 — 1.21 ± 0.08 
to be compared with the global average from 
experimental measurement of /3-decay, 17^4 = 
1.2701 ± 0.0025 [H]. This result provides a 
consistency check of the fitting procedure and 
validates the Bjorken sum rule with an accu- 
racy of about 10%. 

We show preliminary results for the 
NNPDFpoll . set at initial scale Ql = I GeV 2 
together with DSSV08 [5] and BB10 deter- 
minations (Fig. [2]). In general, we can see that 
all PDFs show larger error bands than previ- 
ously estimated, in particular at very small- 
or high-x values, where no DIS data are avail- 
able. This is especially the case for the polarised gluon PDF, which cannot be constrained by 
the available DIS data. We also notice that, at least for the non-singlet triplet, the NNPDFpoll . 
analysis seems to agree better with DSSV08 than with BB10. 

Finally, we compute the first momenta of polarised singlet and gluon PDFs 





A£(Q 2 ) = [ dxAY,(x,Q 2 ) , 
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Ag(Q 2 ) = 


/ dxAg(x, 
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Q 2 ) 




NNPDFpoll.O DSSV08 [5] 


BB10 [2] 


LSS10 [4] 


AAC08 [3] 


A£(Q 2 ) 
Ag{Q 2 ) 


0.31 ± 0.10 0.25 ±0.02 
-0.2 ±1.4 -0.10 ±0.16 


0.19 ±0.08 
0.46 ± 0.43 


0.21 ±0.03 
0.32 ±0.19 


0.24 ±0.07 
0.63 ±0.81 



Table 1: The first momenta of the singlet and gluon polarised PDFs at the scale Q 2 = 4 GeV 2 
in the MS scheme. All uncertainties shown are statistical only. 
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Figure 1: Experimental dataset after kinematic 
cuts for the NNPDFpoll.O analysis. 
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Figure 2: The NNPDFpoll .0 parton set at the initial evolution scale Qq = 1 GeV 2 compared to 
DSSV08 [5] and BB10 [2J determinations. Uncertainties on NNPDFpoll. parton distributions 
are computed at 68% confidence level (see Ref. [5] for details). 

at the scale Q 2 = 4 GeV 2 and compare with the results obtained by other collaborations (Tab.Jlj 
the AAC08 [3] results are given at Q 2 = 1 GeV 2 ). Again, we notice the uncertainties of our 
results: the error on the singlet momentum is between two and four times larger than that from 
other collaborations, while the error on gluon momentum is almost one order of magnitude 
larger. 

More precise determinations of polarised PDFs will have to resort to data coming from other 
processes but DIS, such as open charm and jet production in fixed target experiments or inclu- 
sive jet and W boson production in proton-proton collisions. We plan to extend our analysis 
to these data in the near future, and also to use our PDF set to determine the strong coupling 
constant a s . 
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